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IHTRCD  TACTION 

The  subcritical  assembly  has  been  user?  as  an  educational 
device  to  demonstrate  the  principles  of  nuclear  science  and 

,  and  for  research  in  experimental  reactor 
engineering.  There  are         lvanta~es  in  the  use  of  a 
subcritical  assembly,  Th??  assembly  is  small  in  size  and 
relatively  inexpensive  in  Rati rials,  construction,  ana  re- 
quired l     tus.   Cnly  IK    1  shielding  is  required  due  to 
the  low  activity  involved.   Lattice  p&raaatan  r^n   be  verl 
more  easily  than  in  a  critical  assembly, 

Th   luetic      ions  of  the  assembly  are  time  dependent, 
non-linear  differential  equations  which  may  bo  linearized  in 
the  derivation  of  the  transfer  function.  The  method  of  ex- 
perimental evaluation  of  the  transfer  function  Is  to  vary  cne 
or  more  of  the  parameters  in  such  a  way  that  the  limits  within 
which  the  linearisation  is  v^lid  are  not  exceeded.  This 
variation  of  the  par     rs  is  tailed  "forcing",       foal 

hods  cf  forcing  in  a  subcritical  assembly  are  varying  the 
strength  of  the  neutron  source  or  the  reactivity  in  i  periodic 
manner.       '  I   a  periodic  input  signal  tc  vhich  the  re- 
sponse can  be  measured  is  dc    '  ned  by  the  amplitude  of  the 
random  or  stochastic  variations  in  the  steady  state  flux 
level  due  to  the  inh      statistical  nature  of  the  fission 
prccea  -..   It  is  possible,  in  a  system  v/ith    ■  stochastic 
variations  in  the  neutron  flux,  to  exceed  the  range  of  line- 
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arlty  of  the  kinetic  equations  by  making  the  magnitude  of  the 
input  signal  large  enough  to  obtain  a  readable  response.  A 
technique  has  been       •£,  fron  a     3  of  stochastic 
processes  such  as  Brown!    otion,  tfhioia  -nkes  it  possible 
to  evaluate  th  transfer  function  of  an  assembly.  The 
oehastic  or  random  variatlOB  in  the  flux  is  considered  tc 
the  respon.     the        :c  stochastic  scurce  fcrcing 
due  to  the  statistical  nature  of  the  basic  phenomena  involved, 
This  *     1  of  evaluation  of  the  transfer  function  should 
least  likely  to  invalidate  the  linear  approximations  ma; 


0?  THr  LITFfcATtT 

A  great  amount  of  effort  has  I  in  the  deriva- 

tion and  tv    ticn  cf  crlt.1         ly  transfer  function  . 
One  cf  the  primary  reasons  wss  th«  B#ad  for  an  automatic 
control  system  to  maintain  a  given  power  1  evel  in  a  rnrtr.  . 
The  similarity  of  th-  trir.s:>r  fraction  of  a  reactor  to  th 
Of  a  servomechanism  control  system  had  beas  noted,  and  it  was 
deter  .    that  this  type  cf  system  was  suitable  for  reactor 
control. 

Owens,  Crever,  and     fct  (13 )  in  19^9  proposed  the  con- 
cert cf  tti  ■   |  critical  assembly  as  a  "blaefc  box"  to 
which  a  feedback  control  syste™  eculd  be  applied  in  ordar  to 
maintain  a  desired     r  lavaX*  The  proposal  was  based  on 
the  observation  that  the  roots  of  the  kinetic  equation  of 
assembly  were  negative.  This  indicated  that  the  critical 
assembly  transfer  function  was  th       as  bits  transfer  func- 
tion for  a  minimum  phase  shift  network,  and  that  the  transfer 
characteristic  eculd  be  represented  by  a  Bode  diagram.   In 
that  same  year  Franz  ( [ 0  presented  a  more  detailed  derivation 

the  transfer  funetic  |  that  Hit  standard  tech- 

niques cf  servomechanis     ory  applied.  Harror,      ,  and 
Krucoff  (9)  evaluate      transfer  function  of  the  CP-2  re- 
actor by  sinusoldally  varyiir     reactivity,  and  showed  that 

transfer      ion  could  be  used  in  the  design  of  serve 
loops  or  regulating  ayi      ieh  include  the  reactor. 


ckup  assemblies  her  en  utilized  for  the  determina- 

tion of  the  tra$i*£     >ction  cf   proposed  reactors*      .d, 
Smith,  and  itf.ee  (33    3ured  the  zero  pow       §t*T   function 
of  a  fast  critical  assembly,    -III        :;sembly  of 

trk  III,  b;   determining  the  change  in  power  level  due 
lo  oscillation  of  the  reactivity. 

^inement  in  the  techniques  of  testing  and  the  evalua- 
tion of  transfer  functions  are  constantly  sought,  and  the 
subcritieal  assembly  provides  a  r    :.vely  simple  vehicle 
for  such  investigations.  A  subcritieal  assembly  was  utilized 
by  Axtmann,  D BSSAU  r,  and  Parkinson  (1)  for  reactivity  meas- 
urements and  testing  of  pile  components.  The  transfer  func- 
tion of  the  low-     e  University  subcritieal  assembly  was 
initially  evaluated      cci  (lV)  utilizing  the  method  of 
source  forcing. 

A  new  tech      for  evaluation  cf  the  transfer  function 
of  a  nuclear  reactor  operating  at  a  steady  state  by  an  analysis 

'1c:;,  or  stochastic,  nature  cf  the  neutron  population 
has  been  developed  recently.   cere  (11)  has  shewn  mathemat- 
ically that  the  square  of  the  modulus  of  the  transfer  function 
of  I  reactor  is  proportions!  tc  the  Fourier  transform  of  t! 
autcccrrol  ation  function  which  is  equal  to  the  power  spectral 
density  function  of  the  system.  This  Mil  verified  experi- 
mentally by  Cchn  tk)   and  by  Griffin  m\6   Lundholm  (ft)t     ez 
(lo)  evaluated  the  transfer  function  of  the  Ford  Nuclear 
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letor  and  the  suberit  ly  at     diversity  of 

'  i.^an  by  the  direct  measi  t  o£   th   -^tocorrel ation 

function,  but  the  results  we:-  '    -"Onclusive  in 

the  case  c      a  -ritical  ii  "  y.   MJ  r;  renceg  feu-' 
describing  Hie  evaluation  of  the  transfer  function  by  fteaftt?®* 

the  power  spc  ity  function  pertained  to 
critical  sss'jiblies. 


IONS 

lit  kinetic  equation*  Of  th       I    critical  ftasesl  m 

be  derived  fTOBt  th  ll  filf  KB  eq<  .    (6,   p«    223 

7,   p.   101) 


p 


in  - 1  JljB 


7-0  -  Z90-3   «^|**-§|  (1) 

ch  is  applicabl     that  region  of  a  sufocriti-        ly 

away  from  the  extraneous  neutron  source  where  nest  of  the 

y>e  fission  n outre  r     fee  sy      used  arc  clrfir 

la  th    le  of  symbols,  fbi   first  tera  represents  the  loss 

of  th-  rmal  neutrons  due  to  leakage  from  the  a;    ly;  ti 

eoad  t  r-i  lg  the  less  of  thermal  neutr       to  absorption? 

th«  third  '  Is  th     :ron  production  or  source  term;  and 

t  sic         quation  is  the  non- equilibrium  change 

in  neutron  population  with  time, 

Th     tree  tera  is  mad     of  prompt  Txn.6.   delayed  fission 

neutrons,  and.  nevtrens  rnal  source.   It  has  be- 

shown  (7,  p.  193)  H»t  the  number  of  fast  fission  aetttrei 

produced  is  — -  2.fi  0,  and  further  (6,  p.  226)  that  if  /3 

is  th     ll  fraction  c   I  el  ay    etrons  the       neutron 

source  term  $„   is 
P 

2  t 
Sp  -  (1  »/*)*•  Zj,  0  e~B  (2) 

There  are  six  '       oups  of  delayed  neutrons  that 
account  for  0.61  utrcr     uctien. 


' 
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These  delayed  neutrons  accompany  the  radioactive  decay  of  cer- 
tain fission  products,  called  neutron  precursors,  which  exist 
in  an  excited  state.   If  &    is  the  fraction  of  each  delayed 
neutron  group  then  the  formation  of  the  precursors  of  the 
group  would  be/^  ^"2.0,     When  the  concentration  of  ■  precursor 
group  is  denoted  by  e^  and  th    Mty  constant  as  /L,  the  rate 
of  decay  would  be/l.Cj.  The  net  rate  of  formation  of  pre- 
cursors to  the  1-th  delayed  neutron  group  would  then  be 

4^=  /?nr  za0-  Aici  (3) 

6 
The  total  rttl  of  production  of  delayed  neutrons  is  21  Xe., 

1=1  *  1 
and  when  multiplied  by  the  resonance  escape  probability j  p, 

-B2T 

and  the  slewing  down  non  leakage  probability,  e     ,  the 

result  is  the  delayed  thermal  neutron  source  term 

2      6 

sd  =  p  «rB*"r  H  X  c,  (u) 

a  1=1  l  l 

Table  1.  Table  of  symbols 

B2   Buckling  M .~2 

3         Fraction  of  delayed  neutrons 

Q  ^  Fraction  of  i-th  group  of  delayed  neutrons 

C    Number  cf  delayed  neutron  precursors  in  the 

assembly 
C^   Number  of  the  i-th  group  of  delayed  neutron 

precursors  in  the  assembly 


. 
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1.  CCe   tied) 


c    Concentration  of  delayed  neutron  precursors 
c^   Concentration  of  i-th  group  of  delayed 

neutron  precursors 
B       rmai  neutron  diffusion  ecr  C ''ieiont 
f    frequency 
k      fective  multiplication  factor 

a,  Itipli cation  factor  for  an  infinite  assembly 

L    Neutron  diffusion  length  in  the  moderator 
a  ompt  neutron  li     m   in  a  finite  assembly 

Jj         Pror.pt  neutron  lifetime  in  an  infinite  assembly 
/    Aver     .eutron  lifetime  in  ■  finite  eg    if 
A    Average  decay  constant  of  delayed  neutron 

precursors 
A  *  Decay  constant  of  the  i-th  "roup  of   I  yed 

neutron  precursors 
I      tal  number  of  thermal  neutrons  in  a  finite 

assembly 
n    Number  of  thermal  neutrons  per  unit  vr 
r    He sonar  lity 

0  flux 

3        of  p:     ion  of  "        -.itrons  r 

tal  as  ml   neutron  sou ■ 

contribution  U 


ser. 


gee« 

s<-  e« 


-l 


-i 


..tZ£HSL 
3 


cm< 


r  utrens 

I 

en«  iee« 

neutrcr  3 

3 
j  •  -'  s  ee . 


neutrons 


Table  1.   (Continued) 


^L 

External  thermal  neutrcn  source 
per  unit  volume 

ccntr 

'ibuticn 

neutrrns 
cm.  sec. 

^a 

Macroscopic  absorption  cress  section 

mm     -1 

cm. 

t 

Time 

sec. 

nr 

Fermi  age 

2 

cm. 

V 

Thermal  neutron  mean  velocity 

cm. /sec. 

CO 

Frequency 

radians 
sec. 

The  contribution  of  the  external  thermal  neutron  source  term 
will  be  assumed  tc  be 

S  =    ^L  (5) 

Substitution  of  Equations  2,  U,  and  5  into  Equation  1  yields 

^n   nw2rt   e-  «        ,«         „   \,   <r  rt  -B2T 


■f*  =  D\7«0  -  Ea0  ♦  (1  -/Ok.  2_a0 


e 


6 

-B2r 


(6) 


p  e"    Z_  A4  c.  +   -£- 


1=1 


0  and  Cj.  are  functions  of  space  and  time,  and  it  can  be  sho 
(6,  p.  227)  that  the  space  and  time  variables  are  separable 
allowing  the  partial  differential  equation  to  be  reduced  to 
an  ordinary  differential  equation.   Consideration  of  the 
fundamental  space  mode  only  (6,  p.  227  and  9,  p.  33)  yields 


wn 


- 


1 

the  equality 

V70  *   -  B20  (?) 

'ihe  diffusion  length  is  defined  by 

L2  =  •=£-  (  ) 

*-  a 

Rearrangement  of  the  equation  yields 

D  «  L2Z  (9) 

The  Stan  life  of  thermal  neutrons  in  an  inflnr.         is 
defined  by 

4  s  -V  ■  rV  (10) 

^a 

When  ri arranged  this  tion  becomes 


2.  =  -r~  (ID 

The  flux,  the  effective  M         ion   factor  in  a  finite 
mediu-  ,  and  the         utron  life'    l»  a  firJ 
efined  by 

0  ■  nv  (12) 

,   -B2r 

k  a  -5-5-  (13) 

1  +  IrB 


'*  -    *o 


1    +  ITIT 
ihe  combination  6        ^nations  1.3  ^nd  lU  yields 


■ 


w 
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"  ^fv®27 


w      ~    '„*      bP-r 


Xhe  substitution  of  '.".quations  7,  9,   11,   12,     nd  15  into   I  equa- 
tion 6  yields 


|#  «  -  (L2B2  ♦  1)  *#  ♦   (1   y)k  tp 

4  /  A  (1. 

+  r>  e"'       ZL    A.o.    -     -<z. 
t-1     *   * 

The  combination  c.  -tlon  Ik  with         -  tion  16  and  rearr-  .  - 

is 

2     6 
da  ■     [k(l  -y?)-  l]    f*  ♦  P  e~B  r£   \y     -*.  (17) 

A  i  ~1 

nations  11,  12,  and  15  ftltttttd  i:  le   -nation 

3  the  result  i3 

nk  '  A  -ibr~  -  \h  (lfl) 

P«"J  A 
Equations  1/  and  1  upled  Vrinetic  equations  of  the 

subcritlcal  assembly. 

It  should  be  noted  at  this  point  that  in  arriving  at 

ition  17  th«  as.     ion  has  been  made  that  the  rternal 

source  contributes  a  different  number  of  thermal  monoenr r 

neutrons  to  each  unit  volume  of  the  subcritlcal  assembly, 

the  number  decreasing  with  h      in  the  assembly.  The 

kinetic  equations  can  also  be  deriv    f   considering  tv 

total  neutron  population  ratht     n  the  neutron  population 
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on  a  per  unit   i  >    basis. 

r  that  at  a  particular  time  tfc  bly  contains 
K  neutrons,  and  the  effective  reproduction  constant,  b,  is 
14  than  one.  The  prompt  neutron  lifetime  is  ^  ,  the  pre 
neutron  reproduction  constant  is  ,id  the  rate  of  ill 
of  neutron  pc,  I  ion  in  th-  ly  is  (k  -  1)  M/Jl*%  if 
neutrons  are  not  suppli<  the  bly  from  an  external 
scurce.  Absorption  ■  account  for  tfef 

If  y?  is  i:   .  action  of  i  then  the  pre 

neutron  reprebmetion  constant,  k  ,  can  be  r-     ed  by 
kit   ~/3)>     The  ■  :f  decrease  of  neutron  population  with 

time  is  then    [ires  by 

[Ml  -/?)  -  l]  N//  (19) 

Keutrcns  are  supplied  t  ly  by  an  external  source 

and  six  u.trens.  The  total  number  of 

delayed  neutrens  cf      --th  ^rcup      _1  in  '   'issions 
produf  ibly  at  any  time  is  C,.   If  the  decay 

constant  of  the  nv       ntainin.-  the  i-th  -roup  cf  •' 

trens  is  A*,  the  neutrons  of  that      are  emitted  at 

-   ^i^i        neutrons.  V  fmd   neutrcs 

JtrT 

must  be  corrected  for  fast  i         *,  e     ,  and  resonance 

ipe  probability,  r..  '  the  terms  in  the  non- 

equilibrium  rate  i     bon  for  the  ate 


[ 


|  i  [k(l  »f\*   l]  |*  ♦  pe~B  Tr  \c±+  Se   (20) 

otal  uuabf  &£  produced  j 

unit  ..      *■*§  cf  which  a  fraction  ^f  are  delayed.    .e 
delayed  neutrons  ar«  not  lost  by        cr  absorption  until 
after  I   /  are  emitted  from  the  fission  pre  -.ucts,  and  I 
the  number  of  delayed  neutrons  stored  in  I    'ission  products 

is  gri  than  *jf  by  th*  factor T  g  •  ,j^r  !v:  *  loAlib* 

A  p  e 

rium  equation  for  the  i-th  group  of  d elayed  neutrons  is 

C,  (21) 


u.jsi.  a 


nations  20  and  21  are  the  same  as  liquations  17  and  18  ex- 
y  are         a  the  total  r    p  of  neutrons 
in  the  assembly,  tad  t!  then  accounts 

for  all  of  the  nei    a  supplied  to  the  assembl;  .   he  neutron 
flux     :j   point  in  th       bly  can  be  cons::       o   be 
proportional  to  th*  total  number  of  neutrons  in  the  assembly 
at  a  particular  instant.  Th         ither  sol:  of  e©l 
kinetic  equations  in  the  derivation  of  th*-  transfer  function 
will  yield  the  sa:     ;ult. 


I 


DJ      ION  OF  not  TMNSF 

\o   kinetic  and  the  re     nt  transfer  function 

i  by  varying  the  scurce  or  reactivity  in  a  pre- 
dete:     I  manner.   'Ibis  Variation  is  usually  sinuscid~l  when 
purposely  introduced.   'Jew*.  •/■  2  5  th  scuroe  strength  and  re~ 
activity  nay  vary  even  at  "steady  state"  conditions.    is 
is     U    the  fluctuation  !•  number  of  neutrons  available 

the  assembly  caused  by  the  natural  statistical  fluctua- 
te ns  in  the  rates  of  neutron  absorption  and  fission*  The 
resultant  random  va:i tions  in  th    utrcn  flux      possible 

H  of  th«  transfer  function  by  proper  analysis 
at  steady  s 

ansf er  Function  with  Source  Variation 

Fc       f   state  ccn^i  tiers  (using  zero  subscripts) 
.nations  17  and  il  become 

n       2 
H  »  0  «  [k(l  -*)-  l]  -|  +  pe"B  TL  \c±   *  -<£     (22) 
'       >i         i=l     o    o 


at  B2r   ^*   * 


-b2t  jr    *  i 


o 


(23) 


Consider-  a  small  va  lotion  £^   tc  be  superi-pesed  on  -«£- 

o 

n  ■  n0  +  g  n  (2*0 


ci  =  cic  +  ^ ci 


o 


I     : 


315 

bstitution  of       I    tion  2L   into       I    "Ions  17  *t*&  18  yields* 


(nc  +f n) 


t      6  ^ 

pe"B^T  H  X.Cc,      ♦  £c.)  ♦    -<£     +  ^  (25) 

1=1  o  x 


tndl 


dc<       d  £c4         /?  (n     +  £n)        >  - 

— k  m  Ik  — LA — O -    A    (r     +  oe    )  (26) 

dt         dt  _2T  f  /livcIr       ci;  *w 

pe 

D  the  g  of  Iquatlons  25  and  26  are  expanded,   the 

steac.  onents  of  and  23  ftTt  tOll*€t*4 

and   set  equal   tc  zero,   the  resultant  equations  are 

6 
^  =  ^^n-^-    ^♦^rlyVft        (27) 


-•ansfonsinr   to  Laplace  net  -tlcns  2?  and  28  to 

2       6 
sfa(s)   =  iL&jL   <£"n(s)   -  £4r  fn(s)   ♦  r,e~B  T  21  X^t^i) 

♦  cf^(gi  (29; 


sAe^(s)  »     (*    1^0 att)  -   /\#  d«4<»)  (30) 


£0,(3)   ■  ■■  '  ^     v   on(s)   »    X  &c.  (3) 

Pie 
Rearrangement  of  :  qu.it Ion  30  yields 


f 


A 


~b"t  s  +  xi 


s)  (31) 


',T- 


"2, 


/> 


"i 


T 


ii  A     — _ik 


■ 


S 
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Substitution  of  Equation  31  In  Equation  29  and  colleetl* 

terms  yields 

6 

c5"n(s)  f  s  -  ^4  +  ^4r  -  Z  ^/^  -v  t  f  =  <£^<s)     (3?) 
n\mi    ^    ■    ^    7j^   *_  /*(s+  X1)J  ° 

By  rearrar.      nation  32,  the  transfer  function  becomes 

&$—. ~ ^TTT  (33) 

i=l  s  +  Xi 

By  substitution  of   /  T  for  Z  rr — r~  Equation  33  becomes 

A      i=l  s    Ai 

^44  =  — — *—  (3U) 

^^   s A  k(l  -/)  *  1  -  k  f£~- 

Multiplication  of  the  numerator  and  denominator  by  (s  +  X) 
and  division  by  J*   yields 

I^TTJ     s2  +  (i+ki  +  ^)s^(1  _k)    lw 

The  denominator  is  of  quadratic  form,  and  when  av      values 
of  A  =  0.08,  k  =  0.5,  i*  ■  0.001,  and/^=  0.006U  are 
substituted  in  the  coefficient  of  s,  it  is  seen  that  A  and 
y&   are  negligibl      red  to  — W^  so  that 

<Mg?    s  +  X (^6) 

T^tsT     32  +  (i^)sl  +  jUTjEE 

The  roots  of  the  -re  found  by  use  of  the  quadratic 

formula. 


. 


"\. 


\ 
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s  ■  hn£  f -  i  ± 


2J+ 


1  '  T-rr 


From  the  binomial  expansion  (a  +  b)n  =  a  +  nb where 


1  - 


UA 


1  -  k 


higi.  &£ 


1  -  k 


b   a,  the  radical  expands  to 

By  substitution  of  the  expanded  t     i*or  the  radical,  s 

becomes 


•  •  • «  • 
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First  root 
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When  A  is  considered  negligibl       red  to  the  first  term 
the  resultant  root  is 


s2  =  -  (^ 

The  denominator  is  now  seen  to  be  factorable  to  a  good  ap- 
proximation into 

(s  ♦  A  )(s  ♦  *-*£} 

By  substitution  back  Into  Equation  36  the  result  is 
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(37) 


This  is  the  final  form  of  the  transfer  function.   It  shews 
that  for  a  suberitieal  assembly  one  m*»y  igncre  the  ei     of 


i*v 


'. 


A. 


the  d  •  in  determining  tr  la  Ban 

jtion  tc  a  *  "ion  In  the  sc-rree.     Ih  *mt  ec:  - 

tftl&l&g  tb  rage  delayed  neutron  decay  constant  cancel 

out,   and   thus   the  prompt  neutron  lifet  an  be  replaced  by 

the  average  neutron  lifetime.     The  assumption  that  the  ef 

the  delayed  neutrons  can  be  ignored  is  valid  for  k  #r 
0.995  (1>«     fh*  suberitical  assembly  utilized  in  this  in- 
vestigation has  a  k£,ff  of  .  72,   hence  the  above  assump- 
tion should  be  valid.     The  same  result,   Equation  37  >   is  ob- 
tained f re                ions  1?  and  18  when  I 

containing    ^iC*    IN  deleted  and    the  same  derivation  procedure 
as  above  follow 

Transf  \ction  witfc  -tivity  Variation 

Again  utilizing  1  on  17  and  Equation  18  for  the 

transfer  t  Ion  ar  ~;S  22  ''3  for  the   st 

st-ite  condition,   let  a   ml]    variation     <£  k  be  imposed  upon 

k  so  that 

k  m  k     ♦    £ 
o 

n  =  nc  +    cfn  (3 

where  cfk  is  tha  amplitude  of  the  variation,  and  not  the 
instantaneous  value  of  the  excess  reactivity.  Substitution 
of  Equation  38  into      Ions  17  and  15  yields 
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♦  pe~B  T  Z.  A,(c.      +     £e4)  4     -^  (39) 

1=1     x     *•  x 

dc,        dcfc,  ^    (k  -   <$"k)(n_  ♦  ^n)        . 

t! --ir4, -  <%    °  ,.E2r° Vv^  (U0) 

p/  e 

By  multiplication,    sotting  products  of  snail  numbers   KpSft] 

to  zero,   and  coll cc tin  steady  state  terras  of  Equations 

22  and  23  and   equating  them  to  zoro  one  obtains 

it*  j*        r         r  ~     r 

6 
.  Ilk  ♦  pe-B"rZlL  £%  (Ul) 

d  cf  c,         x?,k_cfn         /2n   Xk 

-rn-1  -       *  °P       ♦       *  c?       -   ^  ,  cfc.  <U2) 

at  *  »B  r  *  -B  T  *      * 

pfe~  p-re"1 

By  transformation  to  Laplace  it  Ion  the   equations  becc 

s^n(s)   =  2s.  ^k(s)  ♦  .£.  cTn(s)   ./~A  <$~n(s)   -/^T£  <£*<•> 

■K  K  X  X 

.  ijiisj.  ♦  pe-E2rZI  X,  cS~c.(s)  CVJ) 

i*  1»1     1       x 


1  /-B2r  f-B2r  *      i 

Solution  of  BqvfttlGfe  Ml  for  cf e^Cs),  substitution  in  Equation 
L3>  and  rearranging  yields  the  tra..sfei  funetloa 
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The  substitution  of  „A/T    for    21   rf — f-     in  Tqu*ticn  U5 

i  ♦>  i=1  s  -    >A 

yields 

A/7 
<£n(s_)        no  X   -   r  +  TTT  (h,s 

vs;       ;c  JL  +  JL.  m  3    4-/3    _  ,A/£. 
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The  multiplication  of  numerator  and  denominator  by  s    ♦    A    and 
collection  of  I  yields 

.£44  =  !k  ,   3(X   -/?)    •»  ^ (1,7) 

o    00  o 

Division  of  the  numerator  ti         tor  by  ^*AQ  yi^l-'s 

O    s2+  (A+__£+£J^)3  ^(—^2) 

This  equation  is  seen  to  have  th     MM  d«no»in«tor  as   Equation 

l-k 
35  which  simplifies  tad  factors   to   (s  +  A)(s  ♦  ■■  •-  ^ ) .      In  the 

numerator  the  footer   (1   ~/3 )   is   seen  tc  be  nearly  equal   to 
one  when  0.0C6L  if   substituted  Tot/3  .  nation  ^8  then  be- 

comes 

£mIm1       V^j 

•    +        y* 

which  is  the  final  form  of  the  transfer  function.  This  result 
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is  the  same  as  Equation  32  w!     the  source  strength  was 
varied  except  ft  tfe   --  itl  ^0/j**     ^he   aver 

cay  const  rt  A  m       icn  of  delayed  neutrons  drop  out 
of  the  final  expression  as  before,  shoving  that  the  transfer 
function  for  a  subcritieal  assembly  is  independent  of  tl 
delayed  n<  "    U  and  dependant  on  the  av      neutron  life- 
tin  i 

odulus  of  the  Transfer  Function 

1-k 


Division  of  the  tr  :  3   r  function,  Equation  ^%  by  — -r 
and  substitution  of  3  *  J  co  ,  yields 


c 


i-k0 

I1  he  res  It  of  normals       uation  '•:  ing  it  equal 

to  T(J  <w)  is 

l**a 

It  is  seen  from  liquation  51  that  the  prompt  neutron  break 

1-k 
frequency,  co    ,  is  equal  to  *y  for  a  subcritieal  assembly. 

co     effectively  deterr?  su 'critical  assembly  if 

yf  is  known,  ane  th       se  is  also  true.  The  m     or 

and  AafU         r<       I  led  by  the  complex  conjugate, 

1  -  IHk^  yitld- 


1  -  ISO* 

"o 

The  modulus  of  the  transfer  function,  or  the  magnitude  of 
Y(j«),  is 

Ml* 


which  reduces  to 


i«J<">  ■-? .  *,  , . a,i  (53) 


■F^y 


The  square  of  the  modulus  of  th«  trans.     unction  |y(J«)| 
is  then 

i"Ko 

The  square  of  the  modulus  of  the  transfer  function  is  t' 
basis  upon  which  tl  s   of  the  transfer-  function  can 

evaluated  as  will  be  shewn  in  the  next  section. 
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The  transfer  function  of  an  assembly  can  be  det    ;ed 
experimentally  by  introducing  sinusoidal  variations  in  the 
input  and  measurin-  th     if  ft&gli   nd  amplitude  attenuation 
with  frequency  of  thv  resulting  sinusoidal  variation  in  the 
output.  Thil  ll  usual 'y  accomplished  by  sinusoidal!  y  varying 
the  stren  th  of  the  neutron  source  or  the  reactivity  and 
noting  the  resulting  variation  in  the  neutron  flux.  Vari- 
ations in  the  neutron  population  due  to  th'  statistical 

UN  of  th  fission  pv#e#fti  cause  a  change  in  reactivity 
and  a  fluctuation  of  the  rower  level  of  a  critical  assembly, 
ihese  rando™  variations  in  the  flux  level  may  be  considered 
as  a  summation  cf  Fourier  s     I  which  transform  Into  Fourier 
integrals  wfc  approaches  Infinity  (2).  Th 

transfer  function  of  an  assembly  could  then  li     ^ly  be  ex- 
tended to  re     it  the  response  of  a  system  tc  a  random  as 
well  as  sinusoidal  input. 

Cons:..   |  voltage  which  is  a  random  or  stochastic  func- 
tion of  til  &   tc  a  small  range  of  frequencies  be- 
tween f^  and  f^.   One  can  dot- [  a  value  which  is  propor- 
tional tc  the  average  power  by  squaring  the  function,  inte- 
grating, and  dividing  the  result  by  the-  time  over  which  the 
integration  was  perform-  .   If  thil  "av  ra:o  ro1  -      :.vided 
by  At  m   f,-  f2,  the  frequency  increment  which  the  average 
power  represents,  the  result  is  th*  average  rower  or  mean 


• : 


square  v       U    olitu"     ft*  unit  frequency.  Thi3  is  rifled 
the  power  speetr  I   nsity  function,  fl    5rm  power  I 
density  function  has  been  extended  to  apply  to  any  function 
which  consists  of  mean  square  signal  awplitv       unit 
frequency  of  the  band  width  lnv     a ted. 

Nuclear  reactions  and  fission  processes  which  give  ft 
tc  the  neutron  flux  in  a  subcritioal  assembly  can  be  con- 
sidered as  a  random  variable  which  has  a  normal  or  Gaussian 

stributicn  about  some  mean  value  cf  frequency  of  occurrence. 
The  power  sp  ctra"     -sity  of  such  a  ran  is  t! 

r«g«  po^rr  (mean  square  signal  amplitude)  rev   unit 
frc  ind  widt       twA    (2).   It  has  the  unie 

property  of  I  i  constant  for  this  I    cf  function  and 

is  defined  as  a  "white  nr.ise".   tfhCB  such  a  random  variable 
is  passed  through  a  linear  systo::.  the  output  is  a  random 
variable  which  has  been  att       by  the  system  character- 
istics (10). 

If  a  system  mefe  ^o  the  subcritical  assembly  is  excited 
by  a  white  nclse,  which  has  a  constant     itudr  it  l! 
frequencies,  one  ran  analyze  the  output  of  tfa       ly  over 
a  range  of  fi         tc  d*1       the  response  or  transf 
function  (8),   It  has  been  shown  mathematically  that  the 
cutput  power  density  spectrum  of  a  chain  reacting  sy-    is 
proportional  to  the  square  cf  the  modulus  of  the  transfer 
function  (11).  Ibis  has  alsc  been  demonstrated  experimentally 


for  Si  rs  (}i,  6 ). 

Figure  1  shews  V   spcrtral  density 

measurement  In  which  it  is  assumed  that  th     otronic 
compener     i        :.  frequency  response  ever  the  range  of 

investigation  and  introduce  no  ftcifl    Under  that 6  condi- 
tions the  syv.  quation  would 

Gyy(W)  .  G^CwjlYCja))!2^:  k|y(j6o)|2    (55) 

re  G  ( oj )  ■  input  eove.r  spectral  density  =  constant. 
G  ( £0  )  =  output  •  cw:  :•  srertral  density. 
Y(j60)  ■  tystca  transfer  function. 
1^0  60)1^  •  tqttarc  of  modulus  of  trarsf     letloa* 
Since  the  frequency  response  of  the  fieaturl     I     nt 
is  net  perfect,  it  alsr  introduer-s  a  transfer  function  into 
the  equation.  This  nay  be  Beg]  erted  ind  the  response  of  the 
measuring  •qnlpa<  at  ew  widened  i:       nt  of  frequency,  if 
the  break  frequency  of  the  equipment  is  beyond  the  range  be:: 
investigated  (S).      fi&al  for*?  of  the  system  equation  was 
ihotm  to  be  ( 

Gyy(w)  =  A  ♦  B|Y(j^>)|2  (56) 

where  A  represents  V  ef  BeamirlB     Ipment  nci 

and  B  represents  the  1        the  assembly  noise.  The  sub- 
stitution of      ion  51'-  into  'quation  56  yi  Ids 
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Figure  I,     Ideal  systora  f©r  determination  of  cutttit  never 

spectral  density  function 
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(6J)  Ri  +  ! — ■—  (57) 

A  sketch  of   I  it ion  5?  is  shewn  in  Figure  2. 
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m  ion  of  ki>  i 

The  graphitt  rated  subcritical  assembly  (iM  utilized 
In  this  investigation  was  constructed  at  lc  -  tat«  University 
in  1957*  Figure  3  is  i  photograph  of  the  suberitieal  assembly 

ch  the  end  cover  r-    d.  The  blocks  were  constructed  from 
AGB  gride  graphite  cylind;         *      :  of  7.0  in.  and 
6  3/q  liW|  and  length  of  6^  in.   The  lower  nine  layers  were 
made  up  of  blceks  with  a  6.0  inch  square  cross  section  and 
rounded  corners  of  radius  3*5  Inches.  The  top  five  layers 
were  constructed  cf  blc     It]    6.0  inch  by  5«0  inch  cross 
section  with  rounded  corners  of  radius  3  3A6  inches.  The 
density  of  the  graphite  was  1.56  gm.  per  cu.  cm.  or  97.3  lb. 

r  cu.  ft.   An  S-J-  inch  lattice  was  formed  by  placing  uranium 
slugs  in  altera- 1   channels  as  shewn  in  Figure  3*  *&•  slugs 
of  uranium  were  1.000  in.  in  diameter  and  8.00  in.  in  length 
and.  encased  in  aluminum  cylinders  0.0^0  In.  thick  and  with 
ends  0.20  In.  thick.  The  outside  dimensions  were  1.080  in. 
in  die         o.^O  in.  in     th.  The  uranium  density  was 
19.0  pu  rer  cu.  cm.  or  11S6  lb.  per  cu.  ft. 

1'he  assembly  was  covered  on  the  top  and  si     ith  a 
0.010  in.  sheet  cf  cadmium  sandwiched  b»  I    B  a  3/8  *n»  sheet 
of  plywood  tttd  t  1/8  in.  sheet  of  masonite.  An  effective 
"black  wall"  for        neutrons  was  produced  by  the  cadmium. 
The  horizontal  distances  bei     the  side  sheets  of  cadmium 
was  60.5  in.  and  62.5  In.  With  the  extrapolation  lengths 
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Fi     ;.  Front  view  of  the  subcritieal  assembly 
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added,  J'  ize  of  tl     sembly  was  approxi    Xf 

62  in.  (ildl  tc  si       l  I  .  (front  to  back). 

IT     (critical  ass  I  3  ^.s 

shown  1   '  ure  3«  ^he  thr  e  *89      housed 

11.5  in.  dee;     -  tanks.  The  water  provi       becticn 
fro      noutrcr         '   **•  mounted  in  the  center  tank 
under  the  pil<  . 

The  neutron  source  consist ed  of  five  individual 
plutcniuE-berylliura  sources  in  contain?  ri  1  in.  in  diameter 

and  1  3/3  in.  in  J    b*      Lad 1^5     sour- 

6 

rated  at  one  curie  emitting  approxin  j  1.63  I  1G     neutrons 

6 
per  sec.     The  tot:  !  itrength  was  B.15  x  10     neutrons 

per  sec.     rihr  five  sources  wer-  on  a  pedestal   in  the 

center  t*Bk  vhlfi  "  r.      The  scurc 

located  on  the  centerli:       I  llreetly  below  the 

bottom  layer  cf  graphite*     Five  holes  W€  1<         the 

fee?  11  '       front  panel   of  thr  assembly  it  one  foot 

tl   intervals.     T1  Lea  V*I      0tt1    to  correspond   to 

vacant  etuu  in  tl      3    ttlee  for  insertion  of  the  BF~ 

neutron  detector*     She  lowest  ■  •  1    ,  cne  foot  above  the  bottom 
of  the  assembly,  will  ed  tc  as  hol<*  ore,   the  hole 

two  feet  in  height  as  hole  two  tad   so  on. 

The  electronic  e:1      I    try  for  Raftenrlng  the  random  varia- 
tions in  the  neutron  flux  is  shown  in  Figure  h.     The  neutron 
probe  used  was  a  wood  Counter  Laboratory  BF-.  detector.     A. 
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Flsrure  *».     Electronic  circuitry 


Nuc"  TJltrascaler  Kcdel  192A  was  used  to  supply 

•h  voltage  tc  the      I  cr  t&d  tc  amplify  the  signal  from 
the  detector.  "12  Direct  Coupl     plifi  r 

With  a  flat  fj    :\cy   response  from  Dr   to  2  nogs cycles  was 
used  tc  fur  J:     plify  the  si       ch  was  fed  Into  a 
Khron-ijite  ultra  low  fi     cy  Band-     Filter,  Kodel  330A* 
The  band  width  of  the  filter  was  ad  '  ■        ble  t     i  a  low 
cut  off  frequency  of  0.02  cps  tc  a  high  cut  off  frequency  of 
2000  cps.   ?•'         'justiGen'    |  accomplished  by  mm  of 
two  dials  calibrated  with  a  logarithmic  scale  reading  direct- 
ly in  cps  from  2  to      A  multiplying    |  e\  tl  of  1/100, 
1/10,  1,  10  and  100.  A  capacitance  f:  ut 

circuit  filtered  out  any  DC  component  of  the  signal.  The 
filter  comprised  an       "y  low  aol   circuit,  with  in- 

loM  than  100  microvolts, 
fill-  1  was  th     d  into  a  two  chane 

word)     k  II  for  furth  r  amplification  and  reccrd- 
•.  Th,  I  cording  pen  sensitivity  was  variable  fzr        .01 
to  10  volts  pox  chart  line,  r:  I  the  chart  speed  was  variable 
from  one  tc  125  mm  per  see.   This  circuitry  was  used  to  in- 
vestigate th*  frequency  range  of  0.02  cps  to  20  cps.   An 
alternate  syr  for  th     yaeney  r     of  20  cps 

tc  100  eps  ftfl  cilitate  data  taking  and 

;sing.    L  of  the  components  after  the  Iftielaar  Chios 
Ultrascalar  were  re?;"'      1th  a  Hewlett-Packard  Uave  analyzer 


. 


'V.  ly  adjusl      "      b  -as 

of  ~s  tc  1V5  epSj  "  "  from 

"  "    :ed  in 

'  jivolts.      Figure   5  is   a        :  equipment 

m  m 

j,  with  the 
probe  on  t'  asse  -y  range 

investir  WM  0.02  cps  to  1  .  try  including 

the  band  p.  '  for   tl  of 

0.02   ^5   t;  •<->.      The  wav.  3  •  -ad   fi" 

tc  100  e  In  o:  introduction  of   spurious 

signals  into  the  system,   th  ':en  fi  -alar 

.hout  fch  !ng  circuit  ictiv l1     *. 

on  fed  into   the  Tektronix  Type  112  amplifier 
input  ai  rlly  with  an 

oscillofjccp I    to   in-  u:  .    .      ickup   sal   no   saturation 

r.     TIm    il       2    '.;as  tl  he 

band  pass  fil  I    a  ©ecu  per  cent  band  pass  such 

.    '     tpf   tc      .  .OS   cps  to  0.16 

0«  output  of  the  bfl  HP  was   recorded 

on  loth   ehftBaa!  -  ush  recorder  for  few©  minutes  for 

each  band   pass  width  setting.     Typical  random  data  for 
representative  bind  pass  li   are  shown  in  e  6. 

Tlv  -   ■    '•■    I  i  alty  ftmcttor  tt  -   th<     »G2  ops  to 

20  cps    I  it  obi  :   for  widtV    I  od 
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ure  ?.   Electronic  equipment 
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ure  6.     Typical  random  data   for  representative  band-pass 

v-idths 
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4f=0.025  cps 
f  =  0.0375  cpt 


I  .-...  ..-.. § 


&f  «0  05   cpt 
f  =0.075cps 


mmmmm  m 


4f  =  0.l0  cps 
f=»O.I5  cp» 


il-020  cps 
f  =  0.30   cpt 


At  =0.40  cpt 
I  '0.60  c  pt 


-U- 


AM  chart    tpeedt    5  mm    par    tec 


41  =  10  tec 


13    r  tc  '                        phy  (12),                        ed 

ordlnates  of  i  pa  trmt                      n         lit  one  s  -- 

intervals.  Oftilf             18   squared,                              ordinates 

and   their  sq"  were   gv                           tly  il                       U        BO 

variance,  (J""",  was  obtained  by  taking  the  diff                          en 

the  average  of  the   squared  vain*                    square  of  the  av 

u.e.      This   is   tqall          t   tc    the   average  cf   I  of 

the  difference  tTw                 ia« 


CT*  *  <y2  -  y2)  =  (y-  y)2 

The  c        ml  lotifttt  Gyy(   )  was  then  c    I  ed 

by  dividing  the  s,nd  r     idth  in  cv^s*     This 

value  cf  r^e^n  squ?r"  a^rliti       unit  free;- 
ey,  or  mean  pe«        ' J-  froqueaey,  for  '.'    nd  width  in* 

ves  ■'     -e? . 

G     (60)   -  O...  3e3n  square  amplitude 

^f.  unit   rrequeney 

e  OOVOl1   spectral  density  for  the  20   eps   to  106  rps 
range  was  obtained  by  sqn  |   the  •  ave 

analyzer  cutrut  aoto?«     The-  meter  road   fliP<  etly  in  HMfl  power 

5  unit  frequency  of  the  band  width  investigated.     By  tqaar- 
illg  the    reading,    the  mean  squai  plltttd*  or  power  per  unit 

frequency  of   th        tad  width  investigated   vol  obtained. 

A  sa:;;l      r-'lerl  ;tion  of   the   type  used   tc    reduce   the  data 
from  the  Brush  recorder  is   included  in  t  ppendix* 


DIJCUI3IC 

The  results  cf  the  data  reduction  are  presented 
graphically  in      3S  7  and  8  with  the  power  sr        nsi- 
ty,  G  (6J),  in  arbitr  ry  units  versus  frequency  in  cycles 
per  second.   hi  data  were  plctt<    ly  ttt        .     ^  to 
10.0  epi  bee*W       ower  spectral  density  a    i ed  to  br- 
corae  ecnstar.t  by  10.0  cps  <md  remained  a  oc     t  frcm  iO»0 
epf  to  IOC  cps.  This  constant  level  r       ts  the  "A" 
in  iquaticn  57   which  is  the  white  noise  of  ■'   Instrumenta- 
tion which  is  tot  primarily  to  the  BF~  detection  chamber. 
The  Inst;    tattoo  ncise  decreas    ' .  th  increasing  height 
In  the  assembly.  This  is  explained  by  the  fact  that  post  of 
the  instru     ticn  noise  1-     tc  the  flux  3     or  count- 
ing rate  in  th      ^tion  chamber  which  i -arrises  with  in- 
creasing height. 

Figure  7  shows  the  change  in  j     I]  atral  density 
function  along  the  axis  of  the  efl    \y,  For  hole  one  (s  = 
1  ft)  the  '       -tral  density  was  determined  to  be  ap- 
prcxir    f   a  constant        at  thl  range  cf  frequence 
investigated.  This  that  the  r      variations  in 

the  neutron  flux  rlly  iue  to  the  flux  frcm  the 

scarce  which  is  "white  ncise".  Hole  two  (z  ■  3  ft)  shows  a 
deviation  from  the  constant  power  spectra;  density  function 
of  the  instrumentation  and  Indicates  an  increasing  contribu- 
tion from  the  fission  flux  and  external  source  which  have 
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Figure  8.  Comparison  of  measured  with  calculated  power 
spectral  density  function 


been  »t1      B  branaf  ft 

'.'.      •    (z  -  3   ft)   is  91. u   a  "      "  •!•» 

le  3    :  '■"        .  i3ity 

'  '  ;    the 

-ion  of  th  red  on  the  flux 

to   a  J  >  r    (2   =  J     ft)  :   a 

may:'  or  attenuation  of  t'  nsity 

by  the  transfer  function  All    -.bly,  five   (z  =  r> 

y  to  approach  a 
•.      Thin  .t  be 

in  tl  -   oi  vtion 

which  r  .3   li  prc- 

root  h«  povmr  or  flux  1-       in 

it at Ion  noise 
■'.  log 

f  of  the  breaks  in  the  eu.   '  s  to        ik  the  in 

of  tfe 

ic   "•   -    la  fro 
data  pT'  in  Fi  7.  "or  helaa   three,   four, 

1th  a  plot  of  the 
•  function  from 

ion  ry  felon  waa  pletis         Lag 

represer  values  of  ^  -       .  .;.        -v-parison 

9  I 


ctral  density  of  h    ."our  to  be  approximately  the  same  &a 
that  for  the  theor-  ig  I  curve. 


CONCLUSIONS 

The  results  were  inconclusiv         the  neutrcn  br« 
frequency  was  not  determiner"  N&t  It  was  not 

rcssibl-  to  eva  lo  ■***  exactly. 

Ihl    r  agreement  between  theory  and  experimental 
results  was  obtain     z  -  k   f t  -   "    t«  At  this  position, 
the  slope  cf  the  attenuated  power  spectral  density  function 

was  approximately  the      li  that  predicted  by  theory.  The 

-k 
value  cf  the  ratio  — y  was  d     ined  to  be  approximately 

.  /0.08  *  5«0.  The  results  Indicate  that  this  ratio  chan 
as  z  :  1  in  value  to  beyond  two  diffusion  lengths  wh 
it  r— tini  fairly  constant  until  masked  by  the  white  noise 
of  th       Ion  chant)--  r. 

Th     lost ion  of  the  transfer  function  of  the  sub- 
critical  assembly  by  analysis  of  the  power  spectral  density 
3hows  promise  and  should  be  pursued  further. 
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8UC  DY 

The  plutorium-beryllium  neutron  sources  utilized  In  this 
investigation  wer.  &n    ■  i  .   "  life  iseaiMy 

and  v-  ipt  for  t  lion  by  the  graphite 

of  till  iy.      In-  by 

encasing  the  sources  in  a  block  of  j  in  ne  -.nd 

surrounc'-:'  -a ted  scurr  ith  graphite  instead  of 

the  water  shield,     ..ere  significant  results  sight  be  obtai. 
by  the  use  of  a  mechanical  or  electronic  device  for  squaring 
and  av  .  of  data  processir.       m  -Id 

allcv;  much  more  rapi  accurate  determination  of  the  pow^r 

spectral  c-nsity,   hence  broadening  the  scope  of  Investigation 
possible  . 
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